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I We propose a method for analyzing Berry phase for a multi-qubit system of superconducting 

O . charge qubits interacting with a microwave field. By suitably choosing the system parameters and 

Q ■ precisely controlling the dynamics, novel connection found between the Berry phase and entanglement 
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^ ; 1 Introduction 

Q^l Berry's phase [1] can be thought of as an adiabatic quantum holonomy restricted to a one- 

dimensional energy eigenspace [2]. The Berry phase has very interesting applications, such as 
CN ■ 

>- ■ the implementation of quantum computation by geometric means [21 IH IS]- Recently, it has 

(N ' 

O '. been recognized that large-scale quantum computers are hard to construct because quantum 

systems easily lose their coherence through interaction with the environment [2]. Researchers 
' have tried to avoid this problem by using geometric phase shifts in the design of quantum gates 

o ; 

0\ I to perform information processing. The robustness of Berry's geometric phase for spin-1/2 



O 



X 



particles in a cyclically varying magnetic field has been tested experimentally [6]. Interesting 
Berry phase results of a composite system have also been presented in [7] and a new formalism 
of the geometric phase for mixed states in the experimental context of quantum interferometry 
has been discussed in Ref. [S]. In another set of experiments, Du et al. [5] performed an NMR 
experiment to measure the geometric phase for mixed quantum states, and Leek et al. [ 10] 
analyzed experimentally the Berry phase for a superconducting qubit affected by parameter 
fluctuations. 

Another important problem associated with quantum computation and information is the 
problem of engineering entanglement in multi-particle systems. This has attracted a great deal 
of recent interest pTl [T2l [T3| [H] . Here questions of how the couplings among the subsystems 
changes the Berry phase of the composite system become of great importance. In this regard, 
the relation between entanglement and the Berry phase has been discussed in solid state systems 
[T5] and in icosahedral Jahn- Teller systems [16]. Most of the earlier works on the geometrical 
phase focus on pure quantum states [TTJ [HI [191 [SHI [21] . These types of systems, however, are 
very unrealistic and almost never occur in practice. In some applications, however, in particular 
geometric fault tolerant quantum computation [5|, mixed state cases are important. From a 



mathematical point of view, Uhlmami was the first to address the issue of mixed state holonomy 
[22| [23] . The difficulty with mixed states is their reduced coherence, which makes any notion 
of a phase more difficult to define and measure. Another challenge for the future is to extend 
this work to more general forms of mixed state multi-qubit systems. 

The aim of the present work is to develop a general treatment of the multi-qubit problem 
and develop realizable procedures to various base spaces. This is dome treating the non- 
degenerate time evolution operator and satisfying parallel transport evolution condition and 
by defining the Berry phase in a new, rather simple, way. In particular, one of the most 
important problems under consideration is how to define the Berry phase of multi-qubit system 
and easy-monitored entangled state with existing experimental techniques [21]. We will show 
that setting up an arrangement in which an entanglement can be measured via the Berry phase 
for superconducting qubits is an interesting task both from the experimental and theoretical 
viewpoints. 



2 Mixed state Berry phase 

The Berry phase has been extensively studied [H [25], [26l [27] and generalized in various ways. 
For example, based on the fact that any mixed state can be represented as a pure one if one 
allows ancillas and optimizes over many purification, the mixed state geometric phase has been 
defined [S], for open systems [2S], and with a quantized field driving [2^ . 

To study the geometric properties of a quantum system [30], we evaluate the Berry phase 
of the system by introducing a phase shift operator 

i?(t) = expH0(t)^V], (1) 

where 0(t) is changed from to 2tt adiabatically. We denote by ip\ip) the creation (annihilation) 
operator of the cavity mode. Then the time independent eigenequation of the system H\xj) = 
'^jlXj), is changed into H'\x'j) = '^'j\x'j), with 

H' = R(t)HRUt) - tRm'^^-J^, (2) 

dt 

and \x'j) = R(t)\Xj)- Hence, the Berry phase generated after the system undergoing an adiabatic 
and cyclic evolution may be calculated as follows 

T 

i,i^) = ^J{x,\R\t)^\x,)dt. (3) 



Whenever a pure quantum state undergoes a parallel transport along a closed path, it gathers 
information on the geometric structure of the Hilbert space in which it lies. On the other hand, 



a state of subsystem is no longer a pure one and to study the Berry phase of the subsystem 
(a non-degenerate density matrix), we have to adopt the definition of geometric phase for a 
mixed state [8]. Our goal in the following is to establish a connection between the Berry phase 
acquired by a multi-qubit system and the Berry phases of the subsystem. 

The geometric phase corresponding to the non-unitary evolution can be defined according 
to the geometric phase of the whole system that evolves unitarily. For example, let us consider 
the following mixed state 

L 

p{t) = j2^kmPkmHt), (4) 

k=l 

where Pk{0) is an orthonormal set of pure states and is the probability, where, p(0) = 
^fc=i ^kPki^) L is the total number of the involved pure states. For non-degenerate C,}. the 
time evolution operator is given by U{t) = Uo{t)^j.e~^'^'' Pfc(O), where au are arbitrary real 
parameters and U^it) is one of the equivalent operators to U{t). When a mixed state given by 
the density matrix p{t) evolves under a unitary operator U{t), the parallel transport evolution 
condition is 

Tr[p,(t)f7(t)f/t(t)] = 0, (A; = 1,2,...,L). (5) 

Based on geometric phase definition for a mixed state [8] and using Eqs. (jl]) and ([5]), for an 
adiabatic cyclic evolution, a new definition of the Berry phase for a mixed states can be written 
as 



7B(r) = argJ^Tr[e,p,(0)f/: 



fc=i 



xexp - / i:i[pMUl{t)Uo{t)]dt 



(6) 



The Berry phase Eq. ([6]) for a mixed state is just an average of the individual Berry phases. The 
Berry phase factor for a mixed state defined by Eq. ([6]) is a weighted sum of the one-particle 
Berry phase factors. In the pure state case Eq. ^ is consistent with the result of Eq. ([3]). 



2.1 Multi-qubit system 

Over the last decade, superconducting qubits have gained substantial interest as an attrac- 
tive option for quantum information processing [311 [32l [33] and Josephson qubits are among 
the most promising devices to implement solid state quantum computation [Ml [35]. A novel 
method for the controlled coupling of two Josephson charge qubits by means of a variable elec- 
trostatic transformer has been proposed in Ref. [36] and the behavior of charge oscillations in 
superconducting Cooper pair boxes weakly interacting with an environment has been discussed 



[38] . Also, the quantum dynamics of a Cooper-pair box with a superconducting loop in the 
presence of a nonclassical microwave field have been investigated in Ref. |39j . 




Figure 1: Schematic picture of the multi-Cooper-pair-box system. The filled circle denotes 
a superconducting island or the Cooper-pair box. It is biased by a voltage Vg^ through the 
gate capacitance Cg^ and coupled to the bulk superconductors by two identical small Josephson 
junctions. The two Josephson junctions have capacitance Cj^ and Josephson energy Eji. The 
total flux is the summation of the static magnetic flux $e, and microwave- field-induced flux 
$/(t) applied via the SQUID loop. 

As an explicit example, we consider a suitable multi-qubit system and envisage a process 
in which the qubits interact with the microwave cavity such that there is no overlap of their 
interactions. Let us start with a short description of a superconducting box formed from a low- 
capacitance Josephson junction of capacitance Cj^ and Josephson energy Ej.. The Josephson 
junction is biased by a voltage source Vg. through a gate capacitance Cg. which is externally 
controlled and used to induce offset charges on the island. The schematic picture of this multi- 
qubit structure is shown in figure 1. The total Hamiltonian of the single-Cooper pair box 
system can be written as [ID] 

where is the creation (annihilation) operator of the cavity mode. In this structure, 

the superconducting island with Cooper-pair charge Q = 2Ne is coupled to a segment of a 
superconducting ring via two Josephson junctions, where e is the electron charge and N is 



the number of Cooper-pairs. We denote by = 0.5(0]^ + the phase difference across the 
junction. The gauge-invariant phase drops 0^ and 02 across the junctions are related to the 
total flux $ through the SQID loop by the constraint 02 — 0i = 27r$/$0; where $o = h/2e is 
the quantum flux. We assume that the structure of the Cooper pair box is characterized by 
two energy scales, i.e., the charging energy and the coupling energy Ej^ of the Josephson 
junction, and we consider that the charging energy with scale Ec = /2{Cg + Cj) is chosen 
to dominate over the Josephson coupling energy Ej and weak quantized radiation field, so 
that only the two low-energy charge states = (|0)) and = 1(|1)) are relevant [H] while 
all other charge states, having a much higher energy, can be ignored. When a nonclassical 
microwave field is applied, the total flux is a quantum variable $ = $e + ^/(t), where $/ is the 
microwave-field-induced flux. If we consider a planar cavity containing the superconducting- 
quantum-interference-device loop of the charge qubit perpendicular to the cavity mirrors, the 
vector potential of the nonclassical microwave field can be written as A{r) = \uk{r)\{Tp^ + V")^? 
where a single-qubit structure is embedded in the microwave cavity with only a single photon 
mode A. Thus, the flux $j can be written as = |$fc|('0^ + ■0), where $fc = <f u^.M. We shift 
the gate voltage Vg (and/or vary $e) to bring the single-qubit system into resonance with k 
photons: E ^ kkujk, k = 1,2,3, .... Note that the charge states are not the eigenstates of the 
Hamiltonian ([8]), so that the Hamiltonian can be diagonalized yielding the following two charge 
states |e) = cost] |1) — sinr] |0) and \g) = sinr] |1) + cos?7 10) with r] = ^ta.n~'^{Ej/2e), where 
e = 2Eci{CgiVgie'^ — (2n + 1)). Employing these eigenstates to represent the qubits, expanding 
the functions cos(7r$/$o) and using the rotating wave approximation, one can derive the total 
Hamiltonian of the system as [311 BT] 



r(i) 



+ E{^^-^Ssin(2,)cos(^) 

m 

^ cos (2r/) i^'g,,^{n)af + H.c] 



fin) a 



Here we denote by a^'* and ai^^ the Pauli matrices in the pseudo-spin basis of the j*^ qubit and 
Qkin), in = Tp^ip) represents the fc— photon- mediated coupling between the charge qubit and 
the microwave field. 

The influence of entanglement on the noncyclic two-particle geometric phase has been stud- 
ied for an entangled spin pair in Refs. [H HJl |33]. It was shown that prior entanglement 
shared between the two spins can change the Berry phase in such an entangled pair. Also, an 
experimental technique for preparing arbitrarily entangled polarization states has been devel- 
oped jS^. The treatment here extends these investigations to explore, in a controlled setting. 



the fundamental features of the relation between Berry-phase and entanglement in multi-qubit 
system. In particular, we consider a process in which multiple superconducting charge qubits 
are placed in parallel and in addition to a static magnetic flux, a microwave field is applied 
through the SQUID loop. Also, each qubit interacts with the microwave field such that there 
is no overlap of their interaction with the field or with each other |35l HHj. In a similar model 
|37], multiple charge qubits are placed in parallel and coupled via a common inductance. 

2.2 Numerical results 
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Figure 2: Berry phase 7b(^) and entropy corresponding to the first instantaneous eigenstate 
versus the dimensionless detuning parameter A/A (A = a/ e^uj/ihCp), for different number of n. 
Here the curves are for n = (solid curve) and n = 10 (dotted curve). The plot was normalized 
and presented in units of vr for the Berry phase. 

In order to understand a number of different cases, we show, in what follows, plots of the 
Berry phase for different values of the involved parameters. As might be expected, the behavior 
of the present system changes depending on the number of Cooper pair boxes (qubits) involved, 
the variation of the system parameters, i.e., the Josephson energy Ej-, gate capacitance Cg- 
and initial state of the field. The Berry's phase defined in equation ([6]) is calculated using 
quantum states evolving in time under the action of Hamiltonian H . If the evolution of only 
one qubit is considered, it is straightforward to write down the formal solution of the problem 
as H\s,g) = fio|s,^), < s < k and = (^^"^ where and f^^"^ are the 

eigenvalues. The entangled states \s,g) and \^±^) are orthonormal and complete. Here ±, 
labels the entangled states which in their bare condition are the ground \g) and the excited |e) 



states, and n, the states of the boson field. The unitary operator U (t) can be written as 

oo 

n=0 /=+,- 
fc-1 

+ ^exp{-imo)\s,g){g,s\. (9) 

s=0 

In Fig. ([2]), we use Eq. IQ and assume that the qubits initially prepared in the excited level 
interact with the cavity under consideration are in the pure number state (Fock state), i.e. 
Pj {0) = \n){n\. In this figure we show the measured Berry phase 'Jb{'^) ^'^d its dependence 
on the scaled detuning A/A, (A = ^/^uj/JhCp), for different value of n. The calculations are 
all carried out with the single photon process and $e = '^'o/2 for all Josephson charge qubits. 
Here two parameters are varied; the number n and the dimensionless detuning parameter A/ A. 
The measured phase is in all cases seen to be exponentially decay with the development of the 
detuning. Detuning is expected to have the same influence on Berry's phase and entanglement. 
This characteristic robustness of Berry phases may be exploitable in the realization of logic 
gates for quantum computation. To compare the Berry phase with entanglement, we use 
the von Neumann entropy, S{t) = Tr[p (t) In p (t)]. The quantum dynamics described by the 
Hamiltonian ([8]) leads to an entanglement between the fleld and the qubits. We can express 
the entropy S{t) in terms of the eigenvalue Tj(t) of the reduced density operator [21], as 
S(t) = — Tj(t) In Tj(t). For a single qubit, disentangled pure joint state leads to S{t) is 
zero, and for maximally entangled states the entropy gives In 2. In our consideration to the 
behavior of the entropy as a function of the detuning parameter A/A. When we take A/A = 0, 
we get almost maximum value for the entropy. As the detuning parameter is increased, the 
entropy as well as the Berry phase are decreased. Further increasing of the detuning leads to 
vanishing of both 7b(t) and S(t), which means a completely pure state is reached (see flgure 
2). 

In the two qubits case, the eigenstates can be written as 

4 

IV'i^'^) = E«?'^HI*^)' (10) 

i=l 

where = |e, e, n), |e, g,n + k), \g, e,n + k), \g, g,n + 2k)), which have an additional phase 
comes from a geometrical feature [H], i.e., the Berry phase is given by 

7„ = 27r(|aiH|2-|a4Hn. (11) 

For the density matrix p{t) = \ipl^^){ipl^^\, which represents the state of a bipartite system, 
concurrence is deflned as 

C{p) = max[0, Ai - As - A3 - A4], (12) 



where the Aj are the non-negative eigenvalues, in decreasing order (Ai > A2 > A3 > A4), of 
the Hermitian matrix T = a/ ^TpP^Tp ^'^^ P = {'^y ® ^ y) P* i^y ® ^y)- Here, p* represents the 
complex conjugate of the density matrix p when it is expressed in a fixed basis and ay represents 
the Pauli matrix in the same basis. The function C{p) ranges from for a separable state to 1 
for a maximum entanglement. Using Eq. (|TT|) . the corresponding concurrences are given by 

Cn = 2max(0, \ai(n)as(n) — a2{n)\'^). (13) 

The concurrences Cn range from (an unentangled product state) to 1 (a maximally entangled 
state). 
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Figure 3: Relations between Berry phase 75(t) and concurrences C„ of the first eigenstate 
varying n from zero for 10 and for different values of the detuning parameter A/ A, where, the 
solid curve indicates the off-resonant case A/ A = 0.3, while the dotted curve shows the resonant 
case A/A = 0. 

Results are reported in Fig. ([3]) for the relationship between the Berry phase and concur- 
rence. From this figure, we confirm numerically that there is a linear relationship between 
Berry phase and concurrence which means that Berry phase is maximum when an eigenstate 
become a maximally entangled state. In a special case, if we consider an entangled Bell state 
given by 

\i^Beu) = ^i\e,e)±\g,g)), (14) 

the Berry phase factor gives exactly the measure of formation of entanglement which is usually 
given by the concurrence [49]. Finally, we would like to make a few remarks on the relation 
between Berry phase and entanglement. Berry's phase has a classical analogue: Hannay's 
angle |50| is a phase effect in a classical periodic system that depends on adiabatically changing 
parameters. When the Hannay angle of a system depends on its action, the corresponding 



quantum system acquires a Berry phase during the same cychc evolution [51]. Also, Aharonov- 
Bohm effects have no classical analogue, but we may treat it as an example of Berry's phase. 
More generally, however, the Aharonov-Bohm and Berry phases can combine in a topological 
phase [52] • 

3 Conclusions 

Based on the above general analysis, the essence of the Berry phase may be summarized as 
follows: The present study obtains explicit results for the Berry phase for a multi-qubit system 
(superconducting charge qubits). We considered the charge-qubit with a SQUID loop and used 
the microwave field to change the flux through the loop. We have illustrated the relation be- 
tween the Berry phase and entanglement by examining different examples in which the Berry 
phase behaves similar to the entanglement. We have shown an interesting phenomenon of de- 
layed Berry phase due to larger detuning that initially entangled junction and field become 
sparable after a finite time. Thus it provides an excellent basis for a further analysis of the 
interplay between Berry phase and entanglement. Finally, one may say that, the supercon- 
ducting qubits model appears quite promising both as a nice theoretical tool, with unusual 
access to exact analytical developments, and in view of physical implementations in quantum 
information processor, including the realization of complex single-qubit manipulation schemes 
and the generation of entangled states. Also, this would help elucidate the relation between 
multipartite entanglement and the Berry phase. 
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